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Abstract 

In this article we present an effective Hamiltonian approach for Discrete Time Quantum Random 
Walk. A form of the Hamiltonian for one dimensional quantum walk has been prescribed, utilizing 
the fact that Hamiltonians are the generators of time translations. Then an attempt has been made 
to generalize the techniques to higher dimensions. We find that the Hamiltonian can be written 
as the sum of a Weyl Hamiltonian and a Dirac comb potential. The time evolution operator 
obtained from this prescribed Hamiltonian is in complete agreement with that of the standard 
approach. But in higher dimension we find that the time evolution operator is additive, instead 
of being multiplicative like that of Ref. [1]. We showed that in case of two-step walk, effectively 
the time evolution operator can have multiplicative form. In case of a square lattice, quantum 
walk has been studied computationally for different coins and the results for both the additive and 
the multiplicative approaches have been compared. Using the Graphene Hamiltonian the walk has 
been studied on a Graphene lattice and we conclude the preference of additive approach over the 
multiplicative one. 


1 Introduction 

Quantum walk (QW) is a formulation implementing quantum principles on the classical random walk 
problem. The classical version in one dimension can be described with a two state coin. As the particle 
reaches a lattice point, the coin is tossed and the particle moves towards either +x or — x direction 
based on the outcome of the toss. But unlike the classical, if we include quantum concept, the particle 
would be in a superposition of +x and —x state in general and the random walk acquires a new nature 
as described in Ref. [2], [3], [4], [5], [6] and others. Quantum walk on a graph with different coins has 
been studied in Ref. [7] , [6] , [8] , [9] in great detail. Quantum random walk has turned out to be very 
useful in developing algorithms for quantum computation as studied in Ref. [10], [11], [12], [13], [14] 
and [15]. 

In Ref. [5] , [4] , the problem of quantum walk on a line was studied in Fourier space and an analytic 
solution was obtained in the asymptotic limit. The effect of absorbing boundaries was discussed in 
Ref. [5],[16]. In case of classical random walk, the distribution attains a Gaussian nature asymptoti¬ 
cally in time with a central peak and the variance is proportional to time; whereas in case of quantum 
walk, the distribution spreads out from the centre and the variance is proportional to the square of 
time as discussed in Ref. [4] , [6] . Quantum walk in two dimension has been discussed in case of square 
and triangular lattice with various coins like Hadamard, Grover, DFT etc. in Ref. [17], [18], [1] and 
others. In this case also the distribution of the probability of finding the particle spreads away from 
the centre in general. 

In recent years, experiments have been performed to implement discrete time quantum random 
walk in case of NMR Quantum Information Processor cf. Ref. [19], [20], trapped atoms cf. Ref. [21], 
ions cf. Ref. [22] , [23] and photons cf. Ref. [24] , [25] . The use of Dirac-like Hamiltonian in the context 
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of unitary quantum cellular automata was discussed in Ref. [26]. Recently the Hamiltonian of the 
Discrete Time Quantum Random Walk (DTQRW) has been studied in Ref. [1] . 

The material in the article is presented in the following way. In Sec. 2 we briefly introduce the 
concept of the quantum walk. In Sec. 3 we give a prescription to formulate the Hamiltonian of a 
QW system from a physical point of view and thereby obtain the time evolution operator from that 
Hamiltonian in Sec. 3.1. We extend our Hamiltonian approach for two dimensional lattice in Sec. 3.2 
and it is found that the evolution operator is additive in nature instead of being multiplicative as 
mentioned in Ref. [1] . The other difference with Ref. [1] is the way in which the effective Hamiltonian 
is perceived. In the above reference the effective Hamiltonian does not have separate information 
about coin operation and lattice translation, where as in our approach the effective Hamiltonian 
shows a clear distinction between the two separate operations. If we take the evolution operator to 
be multiplicative, the walk becomes a two step walk; i.e. the particle must move along Y direction 
after moving along X direction and vice versa. Two consecutive moves along the same direction is 
forbidden. But in our approach, it becomes a single step walk, i.e. the particle is allowed to move 
along the same direction in two consecutive steps. No path is forbidden. This notion of two-step and 
single step walk has been discussed in detail in Sec. 3.2. In Sec. 3.2.2, we compare the results of both 
approach computationally in case of a square lattice using different quantum coins. In the Sec. 4, 
the quantum walk on a Graphene lattice has been studied starting from the Graphene Hamiltonian 
described in Ref. [27]. In this case also we have presented a comparative study of the additive and 
multiplicative approach. 

2 Concept of Quantum Walk 

In case of quantum walk, the particle can be in a superposition of the available basis coin states 
instead of a definite one as in the case of a classical random walk. In one dimension this is equivalent 
to a two level problem. The two basis states of the coin can be represented by the kets j-|-) and j—). 
We assume that j-|-) is a coin state for which the quantum particle moves in the —x direction where 
as for the other coin state |—) the particle moves in the +x direction. A general coin state can be 
represented as 

lx) = c+l+) + c_l-). (1) 

As soon as the particle reaches any lattice point, the state of the coin is changed; e.g. if the coin 
state of the particle was in an eigenstate (1±)), it changes to a superposed state of them. To include 
superposition we need an operator which operates on j-|-) and j—) and results in a superposed state. 
Hadamard coin described in Ref. [4] is a commonly used mixing operator in this regard. The effect of 
Hadamard operator on a coin state can be described as 

Hc\+) = ;^(l+)+ ]-))) 

ftl-> = ^(W-l-)). 

In the following section, we discuss the Hadamard walk in brief in position space formalism. 

2.1 Position Space Formulation 

Following Ref. [5] and [4] , we consider the basis state of the particle as a product state of the chirality 
and the position basis states : j±)jx). Therefore we can define a translation operator T similar to 
that of Ref. [1], as 


( 2 ) 

(3) 


r(n) = l+)(+| (g) jn - l)(nl-h l-)(-| (g) jn-h l)(nl; 

where n is the index of the lattice points. Here we have assumed that at each step the particle moves 
from one lattice point to the nearest lattice points only and the lattice spacing is unity. Now in a 
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basis of 1+) (8) Iat and |—) (8) Iat, T{n) has the form, 


T{n) 


f\n — l)(n| 0 \ 

y 0 |n + l)(n|y ' 


(4) 


Here, Iat ^ is the identity matrix in N dimension whereas N denotes the total number of discrete 
lattice points. Hence the time evolution operator is given by 


W{n) = T{n)Hc 

1 /|n —l)(n| |n —l)(n| \ 
^\^|n+l)(n| -|n+l)(n|y 


Here He is the Hadamard coin. In general, in place of He, we can use any other quantum coin. The 
above one is the expression of the time evolution operator if the particle was at the n-th lattice point. 
Since the particle could be in any of the lattice points, the total time evolution operator will be 


W = ^IT(n). 

n 


( 6 ) 


This is the time evolution operator for 1-D walk in case of position space formulation. 


3 Hamiltonian Method 

In this section we develop the Hamiltonian formulation of quantum walk. In case of quantum walk the 
time evolution operator arises from two processes. First at the lattice points, the chiral state of the 
particle changes and for that, at a given lattice point n, one can write the chirality flipping operator as 
ITi(n) = S ® \n){n\. Here 5 is a 2 x 2 unitary operator. Then one uses a translational operator T{x) 
which moves the particle from one position to another position in space. Clearly, T{x) is a continuous 
operator acting at any point x in space (not necessarily lattice points). In case of the lattice points, 
the time evolution operator for a single step is given by, W{n) = T(n)Wi{n). At places other than 
the lattice points T(x) can itself act as the time evolution operator since at those places, the chiral 
state of the particle remains unchanged. This perception of W{n) immediately leads to the concept 
of an effective Hamiltonian for the system which acts as the generator for time evolution. Like any 
other traditional Hamiltonian this effective Hamiltonian should be Hermitian. 

Since the one dimensional quantum walk problem is a two level problem, we expect the Hamiltonian 
to be proportional to a 2 x 2 matrix. We expect the wave function of the particle to be just translated 
in between two lattice points remaining unchanged in its form. Hence we consider the Hamiltonian 
proportional to kv wherep = hk and v is the effective velocity of the quantum particle, not to be 
confused with the Lagrangian velocity. A straightforward calculation shows that when the Hamiltonian 
kv operates on some wave function, it just translates the state keeping the form of wave function 
invariant. Hence the Hamiltonian representing translation must be proportional to kv. So we expect 
the Hamiltonian to be of the form akv (8) 1 in between two lattice points and (3 (8) + akv (8> 1 on 

the lattice points. Here, \n){n\ is used to make sure that (3 acts only at the lattice points. The term 
akv is valid at all positions in space (even outside the lattice points) and hence we get the term 1. 
Here both a and /3 are 2x2 Hermitian matrices. Now we can write the translational operator T{x) 
(which displaces the particle from a point x to x ± Ax, Ax being the separation between two points 
in space, not necessarily the separation between two lattice points) as 

rrt \ (\^- Ax)(x| 0 \ . 

= ( 0 + 

^In our convention, (where n £ N, the set of natural numbers) denotes n x n identity operator in finite dimensional 
space, 1 denotes infinite dimensional identity operator acting on the space spanned by the coordinates of all the points 
in space and T denotes infinite dimensional identity operator such that X = 1„ (g> 1 
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If we denote the Hamiltonian corresponding to the translation operator T as Ht < 8 ) 1, then following 
Appx. A.l, 


Ht < 8 ) 1 = —kvaz ■ ( 8 ) 

In deriving the previous expression, we have taken h = 1 and have defined v = ^ as the effective 
speed of the quantum particle. Throughout this article, we will always consider h = 1. The derived 
Hamiltonian Ht ® 1 is valid at all positions, whether the position is a lattice point or not. 

If the Hamiltonian corresponding to the coin operation, Wi{n), is denoted by Hs0\n){n\, following 
Appx. A.2, we can write, 

Hs = (9) 

Here, Ar is the time interval to flip the spin, i.e. it is the time interval the particle stays at the lattice 
point n. It is assumed that the particle stays for a very short time-interval at the lattice point n, 
hence At is much smaller compared to the time scale in which the particle is translated from one 
lattice point to the other. We can write the total Hamiltonian as 

H = Ht < 8 ) 1^ (g) |n)(n| . (10) 

n 

Therefore, following Appx. A.3, for any point x the position space representation of the Hamiltonian 
is given by 

H{x) = Ht + '^Hs5{x — n) . (11) 

n 

The Hamiltonian has a continuous term —azkv and a Dirac comb term. This can be interpreted as 
Weyl Hamiltonian with a Dirac comb potential ^ 


3.1 Time evolution operator from a Dirac like Hamiltonian in 1 D 

3.1.1 Formulation and numerical simulation 

Let’s consider the generalized Hamiltonian as derived in Eq.(lO). Therefore, we can write our time 
evolution operator at the n-th lattice point as 

= T{n){S®\n){n\) . ( 12 ) 

We can write down the time evolution operator here as the product of the two terms in this way 
because the term /S.t[Hs,HT\ ® goes to zero due to the following reason. When the particle 

is at some lattice point for the time interval At, Ax ^ 0 and lim^a-^o H^At = 0. Further, when 
the particle is not at a lattice point, H 5 = 0 by our formulation. So we get rid of the commutation 
term. Physically we can say that the particle would stay at a lattice point n for an infinitesimal time 
interval At. In that time interval, the effect of Ht can be neglected. Equivalently, we can say that 
during the chirality flip, the translation is almost zero. Therefore the effect of H^At can be neglected. 
The particle takes some finite time At to reach the next lattice point and when it is in between two 
lattice points, Hs does not have any effect on its state. Hence we can consider the action of the time 
evolution operator for a particular lattice point for a finite time interval At, as an operation to first 
mix the coin states of the particle within an infinitesimal time interval At, followed by a translation 
occurring in a time interval At — At towards a particular direction depending on its chirality. While 
deriving Eq.( 12 ), at the final step we have assumed that Ax = At = 1. 

In general, for a 50(2) matrix S = ( 1 1 we get the spin Hamiltonian, He = 

\smtl cos j 1 ST u 

^At a first glance, Eq.(ll) may look like dimensionally inconsistent but it is not, since while converting a functional 
from its discrete to continuous form we must divide the discrete form with a weight factor of dimension same as that of 
the continuous variable. In this expression we have used a weight factor of unit length. 
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Figure ll in this figure the particle started quantum walk from 
the lattice point 0. After 500 steps probability distribution of the 
outcome is shown in the figure. 


Figure 2) in this figure the particle started quantum walk from 
the lattice point 0. After 500 steps probability distribution of the 
outcome is shown in the figure. 


For a SU{2) matrix S = 


cos 0 isin0' 
, i sin 6 cos 9 , 


, we have Hs = 


At 


previous expression. 


then we get, ^ 


1 


and thereby get Hs = —j^crx- The result of 


1-D quantum walk for this coin is shown in Fig. 2. If we choose the Hadamard operator as S' = He = 
^.„ ... . „ ^ ^ ^ 


1 


then the spin Hamiltonian Hs = 


72 \ I -1 I ’ - 2AT I J_ -1 - J- 

\ / \ 72/ 

walk with the Hadamard coin is shown in Fig.l. For both Fig.l and Fig. 2, we have chosen our initial 

chiral state to be |y;) = -^(1+) + |~))- The particle starts at the middle position of the lattice points, 
so its initial position state is |^) = |0). Therefore, the total initial state of the particle \7) = |x) ® 10- 
The results of Fig. 1 and Fig. 2 are the standard results of one dimensional quantum walk well studied 
in the previous literatures. In rest of the cases of this article, we will generally choose that the particle 
starts from the middle point of the lattice structure unless otherwise stated. So in general, only the 
initial chiral states will be mentioned. Always the initial state of the particle is chosen to be a direct 
product state of the coin state and the position state. 


The result of quantum 


3.1.2 Calculation of variance 

We calculated the variance of this 1-D walk for different no. of steps numerically. The result shown in 
the Fig. 3 is for the Y coin and an initial chiral state |x) = ■^(l+) + l~))- The plot can be approximated 
as a polynomial of the form 0.47175578 -|- 0.00091395f -|- 0.29289026t^. Clearly the leading order term 
in the variance unlike the classical walk as already pointed out in Ref. [6] and others. 
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Figure 3: Plot of variance vs. time in case of 1 D quantum walk. The errorbar is with respect to the 
quadratic polynomial is multiplied by 10^. The Y coin has been used to generate the plot. 

3.2 Two Dimensional case for a square lattice 

In Ref. [1], while formulating random walk on a square lattice, the Hamiltonian has been taken in 
additive form and the total time evolution operator has been taken as the product of the evolution 
operator acting along each axis. It implies that the particle has to move in x — y or y — x order only. It 
does not have the choice to move along the same axis in two consecutive steps. This approach can be 
thought of taking two steps (in x — y or y — x order only) in each operation of the evolution operator. 
So we call it two step approach in the following sections of this article. But in our approach, the 
particle is allowed to move along any axis in each step. Hence in this single step approach there is no 
forbidden path for the walk. 

3.2.1 Formulation 

In random walk, a particle can move along any of the paths available to it from a lattice point. So in 
case of quantum walk the dimensionality of the coin (i.e. the dimension of the Hilbert space on which 
the coin acts) is determined by the number of paths the particle can choose. In all dimensions higher 
than one, we have two types of chirality of the particle for the approach described in this article. One 
is to choose which axis the particle would move along in the next step and the other to decide whether 
the particle would go to the positive or to the negative direction of the chosen axis. 

In case of a square lattice, the particle has two choice of axes at each lattice point. In this case, 
analogous to Eq. (11), we can say that the Hamiltonian at a point (x,y) for 2-D walk will have the 
form 


H (x, y) = ^ Hs6{x — n)6{y — m) + Ht <8) 1 . 

n,m 


Here, 


Hj' 


Hs = In {Sx ® Mi + Sy® M 2 ) 



+ kyV, 


yUy 



Where Sx and Sy are given as. 




(Sll 5i2\ 


Sy = 




(13) 


(14) 

(15) 


(16) 


and 
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We define S = Sx + Sy as the chirality operator for choosing the X or the Y axis. Mi and M 2 
are chirality operators for choosing positive or negative direction of a particular axis. Now, if we 
choose discrete set of lattice points such that the position eigenstate is |n,m), with the condition 
{n,m\p,q) = dnpSmq {n,m,p,q are integers), then following similar approach as shown in Sec. 3.1.1, 
we can write (Appx. B), 

Tnm {{Sx ® Ml + S'y (g) M 2 ) (g) |n,m)(n,m|) , (17) 

0 0 0 \ 

\n + l,m){n, m| 0 0 

0 |n, m — l)(n, m\ 0 

0 0 |n, m + l)(n, m\J 

is the translation operator. The result of Eq.(17) shows that the time evolution operator is additive in 

nature. This clearly differs from the two-step ^ approach of Ref. [1] where it is assumed that the two 

dimensional time evolution operator is the product of one dimensional time evolution operator in the 
X and Y direction. Now for the total chirality operator M = Sj, (g) Mi + Sy<Si M 2 we can use different 
types of mixing operators, e.g. Grover, DFT or the higher dimensional Hadamard coin as discussed 
in the next section. 


IRnm — 


where 


T — 

-L n.m. — 


/|n — 1, m)(n, m\ 
0 
0 

V 0 


3.2.2 Plots for different coins 

Case 1 : First let’s choose the coin as the direct product of two Hadamard coins. 



Figure 4: Probability of finding the particle after 
250 steps starting from the origin. This plot is with 
the two-step approach using the coin of Fq.(18). 


Figure 5: Probability of finding the particle after 
250 steps starting from the origin. This plot is with 
the one-step approach using the coin of Eq.(18) 


In that case, 

/l 1 1 

M=-; \ 

2 11-1 

_ VI -1 -1 

® In the two-step approach if the particle moves along x axis at one step, the next step must be along y-axis. So we 
can interpret the method as follows: first Wx operates on the particle and displaces it along ^-direction. When it reaches 
to the next lattice point, Wy operates. So in this case the time evolution operator is different at different lattice points. 
The total time evolution operator for the lattice is the sum of the time evolution operator at all the lattice points. If the 
particle starts to move along X direction from the point (0,0), at all the even points (points for which x -\-y \s even) the 
operator Wx will act while at all the odd points (points for which x + y is odd) the operator Wy will act. Although the 
operator is not actually multiplicative, the net effect in two-step can be described by WyWx because after evolution by 
Wx the particle would reach next lattice point where the evolution is governed by Wy. This can happen if at the points 
where Wx is operating, kyVy = 0 and at the points where Wy is operating, kxVx = 0. In this way, we can reproduce the 
two-step approach of Ref. [1] starting from our Hamiltonian formulation of single-step walk. 


-1 

-1 

1 / 


(18) 
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We choose our initial chiral state to be |x) = (8) -^(1+) + i|—)) . From now onwards 

|1), |2) will denote the chiral state of the particle to choose the X and the Y axis respectively and 
|+),|—) denote whether the particle will move along the positive or negative side of a chosen axis 
respectively. The result of the walk governed by this coin is shown in Fig. 4 and Fig. 5. The spin 
Hamiltonian for this coin is given as 


Hs = 


Case 2 : 



1 

-1 

1 

1 

1\ 

TT 


1 

-3 

1 - 

-1 

4 At 


1 

1 

-3 - 

-1 


[ 

1 

-1 

-1 - 

-1/ 

’ coin 

described in Ref. 

[6] 



f-1 

1 

1 

1 

1 


1 

-1 

1 

1 


11-11 
VI 1 1 -V 


(19) 


( 20 ) 


I 


Prdbabillty 



I 

Via. 


0.405 


0.304 


Probability 


0.203 


250250. 



-250. 


125. ^^^^125. 
250250. 


Figure 6: Probability of finding the particle after Figure 7: Probability of finding the particle after 
250 steps starting from the origin. This plot is 250 steps starting from the origin. This plot is 
with the two-step approach using the Grover coin ^ith the one-step approach using the Grover coin 
ofEq.(20). ofEq.(20). 

and we choose our initial chiral state as |x) = ;^(|1) ~ |2)) ® -^(l-l-) — |—)). The result is as shown 
in Fig. 6 and Fig. 7. The spin Hamiltonian Hs for this coin is given by, 


Hs = 


IT 


4At 


f-3 1 1 1 \ 

I - 311 

II - 31 

VI 1 1-3/ 


( 21 ) 


Case 3 : Now we choose the DFT coin described in Ref. [6]. Therefore, 


fl 1 1 

1 i -1 
1 -1 1 
VI -i -1 


1 \ 

—i 

-1 

* 


{ where, i = V-X ), 


( 22 ) 
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Figure 8 : Probability of finding the particle after 
250 steps starting from the origin. This plot is 
with the two-step approach using the DFT coin of 
Eq.(22). 


Figure 9: Probability of finding the particle after 
250 steps starting from the origin. This plot is 
with the one-step approach using the DFT coin of 
Eq.(22). 


and we choose our initial chiral state |x) = ~ l^)) <8> '^(1+) “ !“))• The result is as shown 

in Fig . 8 and 9. The spin Hamiltonian for this coin is given by 




1 

1 

1 \ 

TT 

1 

-2 -1 

0 

4At 

1 

-1 -1 

-1 



0 

-1 

- 2 / 


(23) 


The plots of Fig. 4, 6 and 8 are the standard two-step quantum walks well studied before in the 
other literatures but they have been shown here to compare the results with that of the single-step 
approach of the corresponding coins. 


3.2.3 Calculation of variance 



Figure 10: Plot of variance with time in case of a square lattice. The red curve shows variance in case 
of the single step approach and the blue one for two step approach. The errorbar with respect to the 
quadratic polynomial in both cases is multiplied by a factor of 10^. The DFT coin of Eq.(22) has been 
used to generate the plot. 
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We calculated the variance of the probability distribution for the square lattice case and result has 
almost similar form with all types of coins we have used. So we mention here only the result obtained 
with the DFT coin. The above figure is obtained with an initial chiral state, |x) = + |2)) ( 8 ) 

■^(1+) + |—)). The red color plot in Fig.lO is for single step approach and can be approximated with 
an approximated polynomial of the form 0.42904373 + 0.08807425t + 0.07993539t^ and the blue color 
plot is for two step approach and can be described with an approximated polynomial 1.09843965 + 
0.17597539t + 0.15954993t^. So we see that the variance of the probability distribution in QRW has 
quadratic nature with time in case of 2D square lattice, similar to the form for ID case as mentioned 
in Ref. [4], [ 6 ]. 


4 Quantum walk on Graphene Lattice 

The Graphene Hamiltonian according to Ref. [27] is given hy H = vpa.p which is valid only at 
positions very close to a lattice point. Taking hvp = 1, we can write H = a.k = Uxkx + (^yky. We 
will see that this Hamiltonian can generate an infinitesimal time evolution operator close to a lattice 
point. To generate a finite time evolution operator we need the Hamiltonian for translation operator 
to act at positions in between two lattice points. 


4.1 Formulation of the problem 

To describe random walk problem on Graphene lattice, we have to consider the three available paths 
at each point. So we need three axes as defined by ni, n 2 and ns as shown in the Fig. 11. We denote the 
directions along these three axes by 1, 2 and 3 respectively. Glearly these three quantities ni, 712,^3 
are not linearly independent but they are dependent by the relations 


X = 


V^(n2 - ns) 


and y = —ni + 


(n2 + ns) 


(24) 


By choosing these three axes we can easily extend our previous approach for square lattice in this case. 
In this approach, the position ket at each lattice point (x,y) is denoted by |ni,n 2 ,n 3 ). For Graphene 



Figure 11: Here Graphene lattice is presented with the three axes ni, n 2 and ns. Some of the lattice 
points are specified on the lattice. 

lattice, we can classify the lattice points into two types. The points for which the sum of ni, n 2 and 
ns is even, are called even points and those having an odd sum of lattice indices are classified as odd 
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points. Let us consider the two nearby lattice points, say the even point (0,0, 0) denoted by point 1 
and the odd point (0,1,0) denoted by point 2. From the point 1, the particle can move only to the 
positive direction of all the three axes (see Fig. 11). But from the point 2, the particle moves only along 
the negative direction of all the three axes. Clearly, the geometry itself constraints the direction of 
movement here. We can describe this by using two kets as |+) and |—) which denotes the chirality of 
going along the negative or positive direction respectively as in the earlier cases. The only difference 
in this case is that the particle cannot be at superposition of |+) and |—), it must be at either of 
them. It is clear from the geometry that if the particle is in |±) state at some step, at the next step it 
must be in |=f) state. For the lattice as described in Fig. 11, for even points, the particle is always in 
|—) state and for the odd points it is always in |+) state as per our convention. Clearly the operator 

^0 


needed to describe this kind of motion is given by 


b 0 


Again, we need a 3 x 3 matrix to make the 


superposition of three chiral states which determine the choice one of the three axes. After choosing 
the chirality, the state must translate. Following Eq.(17), we can write the general time evolution 
operator as VF(ni,n 2 ,n 3 ) = T{ni,n 2 ,n 3 ) {{Si 0 Mi + S 2 0 M 2 + S 3 ^ M 3 ) (g) |ni, n 2 , n 3 )(ni, n 2 , nsl). 
Here ^i, S 2 and S 3 are 3x3 matrices which have only the 1^*, and 3^^^ rows non-zero respectively. 
We can choose any three dimensional unitary operator as S (where, S = Si -|- ^2 -|- ^ 3 ). We would use 
the Graphene Hamiltonian to find out Mi, M 2 and M 3 . We can express kx, ky the components of the 
wave vector k, in terms of ki, k 2 and k 3 in the following way. 


\/3 1 

kx = — {k 2 - k 3 ) and ky = -ki + -{k 2 + ^ 3 ) 


So close to a lattice point, the Graphene Hamiltonian can be written as 


Hxy = {kxcrx + kyay)®\x,y){x,y\ 

_ 0 ikA ( 0 iuP‘k2\ ( 0 iwA:3\\ 

ifei 0 j \^—iujk2 0 j \—ioj'^k3 0 j j 


x,y){x,y\ . 


(25) 


(26) 


where, oj = i = Since the above form of the Hamiltonian is valid close to the lattice 

point, i.e. within a very small region around a lattice point, it would represent infinitesimal time 
evolution of the quantum particle near the lattice point. This Hamiltonian shifts the chirality of the 
particle and translates it infinitesimally. The Hamiltonian for a position away from the lattice point 
is given by Ht, the translational Hamiltonian. So first the Graphene Hamiltonian acts and leave the 
particle slightly away from the lattice point. After that the translation operator takes it to the next 
lattice point. Now following Eq.(26) and neglecting the subscript in Hxy, we can write. 


-iH = -i 


0 iki\ I 0 ioj‘^k 2 \ 

—iki 0 j \ —iu)k 2 0 / I —ioJ^k 3 0 


0 iujk 3 


\ni,n 2 ,n 3 ){ni,n 2 ,n 3 \ 




0 1 -|- iki 

1 — iki 0 


0 


- u:‘^{l + ik 2 )\ ( 0 

\u}{l-ik 2 ) 0 ) \ui‘^{l-ik 3 ) 


uj{l + iks) 
0 


. f 0 1 + OJ + 

+Ml+^ + a;2 


0 


\ni,n 2 ,n 3 ){ni,n 2 ,n 3 \ . 


X — iH/\t = 


—i 


+ 1 ,2 


0 l-|-zfciAxi\ / 0 oj‘^{l + ik2A.X2)\ 

1 — ikiAxi 0 y yoj{l — ik 2 Ax 2 ) 0 y 

0 a;(l-I-ifc3Ax3)\ 

oj'^{l — ik 3 Ax 3 ) 0 y 


\ni,n 2 ,n 3 ){ni,n 2 ,n 3 \ +1 . 
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Neglecting the identity operator above (which is not involved in any physical operations) we get three 
infinitesimal time evolution operators as 

,n 2 ,n 3 ) 

,n2,n3) 

W3{ni,n2,n3) 

In the previous equations we have assumed that the —i can be absorbed in the normalization of the 
initial chiral state of the particle. More over the Axj’s which appear are assumed to be very small 
displacements around the lattice points, much smaller than the lattice spacings. As the particle is 
Axi distance away from a typical lattice point the Hamiltonian is governed by the term Ht, which 
is obtained from the translational operator T'(ni ± 1 , ni ± Axi; n2 ± 1 , n2 ± Ax2; n.3 ± 1 , na ± Axs), 
responsible for translating the particle from the very neighbourhood of a lattice point to the adjacent 
one. When T' operates on (Si (8> Wi + ^2 ® W2 + S3® W3), we can write the total time evolution 
operator in terms of the of an effective translation operator T (77,1,77,2,^3) as 

W (771,77,2,773) = T (ni, 772 ,773) ((<51 0 Ml + S2 0 M2 + S3 0 M3) 0 |77i,772,773)(77i, 772,?73|) ( 27 ) 




0 1 771 - Axi, 772, 773) (771,772,773 1\ 

J?7l + Axi, 772,?73)(77i,772,773| 0 j ’ 

0 a;^|77l,772 - Ax2,?73)(77i, 772,773| 

a;|77l, 772 + Ax2,?73)(77i, 772,773] 0 


0 

07^1771,772,773 + Axs) (77l, 772, 773 ] 


a;|77l,772,773 - AX 3 ) (77l, 772,7731 


0 


where. 



M 2 



and M3 



and T (771,772,773) is the effective translational operator defined as. 


r(77i,772,773) = diag (|77i - 1 ,772, TTa) , |77l + 1, 772 , 773) , 1 77l , 772 - 1, 773) , |77i, 772 + 1, 773) , 1 77i, 772, 773 - 1) , 

|77i,772,773 + 1)) (771,772,773] . 

Here the effective translation operator T is a combination of two translations 


4.2 Plots for different coins 

Now for S, we can choose different types of mixing operators. 


Case 1 : Let’s take S to be the three dimensional DFT coin. Therefore, 


S 




LO = 


-1 + \/37 

2 


( 29 ) 


We choose our initial chiral state to be ]x) = ( 1 ^) + 1 ^) + 1 ^)) ® !“)• results are shown in 

the following figures. The plot of Fig. 12 is generated by taking the additive time evolution operator 
(i. e. single step approach) and for that of 13 , we take the time evolution operator in a multiplicative 
form (which is equivalent to a three step approach here). 

In case of three step approach, we see that the probability of finding the particle collapses to a 
single point only. This is expected from the formulation of QRW, since for the multiplicative form of 

■^The translational operator is generated by two translational operators as follows: 

T(ni, 772, na) = T' (ni ± l,ni ± Aa;i; n 2 ± 1, 712 ± Aa; 2 ; 713 ± 1 ,773 ± A^s) (28) 

xT''{ni ± Axi, 77i; 772 ± Ax2, 772; 773 ± Ax 3 , 773) 

In the above expression, the T” term comes from the Graphene Hamiltonian itself and the T' term comes due to putting 
the extra term involving Ht into the Hamiltonian. 
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Figure 13: Plot of the probability of finding the 
Figure 12: Plot of the probability of finding the par- particle after 200 time steps using the three step 
tide after 200 time steps using the single step ap- approach and the three dimensional DFT coin of 
proach and three dimensional DFT coin of Eq.(29). Eq.(29). The values in the color-bar are multiplied 

The values in the color-bar are multiplied by 50000. i^y 200. 

the evolution operator a particle starting from origin always reaches the point (1,-1,1) for every odd 
step and returns to the origin for every even step. 


Case 2 : Now let’s choose S to be the three dimensional Grover coin. Therefore, 

I” 

179. 



Figure 14: Plot of the probability of finding the particle after 200 time steps using the single step 
approach and three dimensional Grover coin of Eq. (30) . The values in the color-bar are multiplied by 


5000. 


S = 


1 

3 



2 

-1 

2 



(30) 


The result for single step approach is shown in Fig. 14. Here we have chosen our initial chiral 
state to be |x) ~ (1^) + *1^) “ *1^)) ® !“)• three step approach is same as Fig. 13. 

Therefore, we omit it for the sake of brevity. 
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4.3 


Calculation of variance 



Time 


Figure 15: Plot of variance with time for the quantum walk on Graphene lattice. The blue line shows 
the variance for the single step approach and the red one for the three step approach. In case of 
the single step approach, the errorbar w.r.t. the quadratic polynomial is multiplied by 50. Three 
dimensional DFT coin of Eq.(29) has been used. 

We calculated the variance for different time steps for both single step and three step approach 
using the coins of Eq.(29) and Eq.(30). Since the results were similar, we include here only the result 
obtained with the three dimensional DET coin. The plot with blue line in Eig. 15 is for evolution 
operator in the additive form and that of the red line in the same figure corresponds to the time 
evolution operator in multiplicative form. Eor the single step approach, the relation of variance with 
time has an approximated polynomial form 0.34339904+0.03773854t+0.02227997t^. But for the three 
step approach, we see that the variance is zero, which is expected since the particle always reaches a 
single point in this case as mentioned previously. 

5 Discussion 

Our approach to DTQRW sets constraints to the form of the Hamiltonian of the system. We found that 
the QW Hamiltonian can be described by the Weyl Hamiltonian with a Dirac comb type potential. 
In Ref. [1], though it is suggested that the Hamiltonian for QW in one, two and three dimension 
resembles a two component Dirac like Hamiltonian, our form of the effective Hamiltonian is distinctly 
different from the form of the Hamiltonian prescribed in the above reference. The other important 
difference with Ref. [1] for multidimensional lattices, is the way in which the evolution operator is 
perceived. In the referred work these evolution operators are multiplicative while in our approach the 
evolution operator is additive for dimensions higher than one. We have shown that the multiplicative 
evolution operator can be produced from our formulation also in case of two step approach by setting 
/ca; = 0 or /cy = 0 at alternative lattice points. 

The Eourier space formulation of QW as studied in Ref. [4] gives an analytic asymptotic form of 
the wave function of the particle. Hence this method results in a good understanding of the problem. 
In this approach, the method to solve the problem reduces to diagonalization followed by an inverse 
Eourier transformation. It works very well in case of ID and 2D square lattice problem, but in case 
of more complicated lattice like Graphene lattice, this method becomes too complicated. In this 
situation, position space formulation turns out to be more useful. Position space formulation gives a 
formal way for numerical simulation of the problem of DTQRW for any kind of lattice. 

In this article, we have calculated the variance of the probability distribution and found its 
quadratic dependence with number of steps. With a little modification of the standard Graphene 
Hamiltonian, we obtained the quantum walk on a Graphene lattice. Since the Graphene lattice is 
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more constrained, we expect more centred probability distribution which is in agreement with our nu¬ 
merical computation of the walk as shown in Fig. 12 and Fig. 14. Due to the geometry of the Graphene 
lattice, |-|-) and |—) state can not exist for the same lattice point. So there is no interference between 
those states. Still in this case we find the quadratic dependence of variance with time Fig. 15. The 
possible reason for this quadratic dependence may be the interference between |1), |2) and |3) states. 
Use of multiplicative time evolution operator in this case leads to the confinement of the particle at 
a single point and thus no spreading occurs as seen in Fig. 13. Hence we conclude that, the additive 
time evolution operator is preferable over the multiplicative one. 

In short, we obtained an effective Hamiltonian for QW from the evolution algorithms already exist¬ 
ing in the literature. The form of the Hamiltonian became important when the analysis on Graphene 
lattice was presented, because a standard Hamiltonian of similar kind (linear in k) was already ex¬ 
isting there. Using the similarity of the Graphene Hamiltonian and the effective Hamiltonians for 
QW a new algorithm of QW on Graphene lattice was presented. In a nutshell, we have presented a 
new and interesting way of attaining an effective Hamiltonian for QW which reproduces most of the 
results of earlier work in the same direction in one-dimension but differs from earlier results in higher 
dimensions. 
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A.2 Derivation of Hs 

If we denote the Hamiltonian governing the chirality flip as J2m ® |m)(m|, we can write 

e-iErr.^sAr^\m){m\ ^ S <S> H {m\ 

m 

m 

^ (g) |f7,^ = s ® |n) 

^Hs = ^ln5. 

Ar 

Now Hs is Hermitian by construction, since every unitary operator on a Hilbert space can be written 
as [/ = for some Hermitian A. 


A.3 Derivation of H(x) 

The total Hamiltonian H is given by, 

H = \n){n\ + Ht <8* 1 

n 

Therefore, 

H{x) = (I 2 ® (a:|) H (I 2 <81 |a)) (|a) is some ket G IR, the position space.) 

= HsS{x — n){n\a) + HT{x\a) 

n 

= (^Hs5{x-n) + Ht^ {x\a) 

^H{x) = ^ Hs5{x — n) + Ht ■ 

n 

Here in the 2'^^ step we have used a proper unit weighting factor. 
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Appendix B Derivation of Evolution operator from Hamiltonian in 
2D square lattice 


We have, 


H{n,m) = Ht 0 1 + Hs 0 \n,m){n,m\ 

Now following the logic as described in Sec. 3.1.1, we can say that [Hs, Ht] At = 0. Now, 


Hs = —In Ml+ Sy0 M 2 ) , 
At 


^—iHsAT<S>\n,m) {n,m| 


Again, 


g "“a—'( 1 ^ (g) 777,^^72^ = (^Sx 0 Ml + Sy 0 M 2 ) 0 \n, m) {u, m\ 


Ht = -IkxVxT^ 


Therefore, 

^T{x,y) = 

Therefore, we have. 


0 OJ '^ \ 0 e 

AkxAx 


1 0\ _ 0 \ 0^ 0 

—ikxAx I 


,0 1 


^gifc^Ax Q Q Q 

0 0 0 
0 0 0 
V 0 0 0 } 


^ 0 e~^>^y‘^y ^ 

® \x,y){x,y\ ■ 


® \x,y){x,y\ 


W{n,m) = T{n,m) {{Sx'S i Ml + Sy'S> M 2 )'Si \n,m){n,m\) 
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